Rogue waves emerging from the resonant interaction of three waves 
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We introduce a novel family of analytic solutions of the three- wave resonant interaction equations 
to the purpose of modeling unique events, i.e. "amplitude peaks" which are isolated in space and 
time. The description of these solutions is likely to be a crucial step in the understanding and 
forecasting of rogue-waves in a variety of multi-component wave dynamics, from oceanography to 
optics, from plasma physics to acoustics. 

PACS numbers: 05.45.Yv, 02.30.Ik, 42.65.Tg 



Introduction. The nature of rogue waves, mostly 
known as oceanic phenomena responsible for a large num- 
ber of maritime disasters, has been discussed in the liter- 
ature for decades [ll-Q . A number of various approaches 
have been suggested to explain the high-impact power of 
these "monsters of the deep" [5], which appear visibly 
from nowhere and disappear without a trace. Theories 
may differ depending on the physical conditions where 
these waves appear [6|, 0] . 

As a matter of facts, a comprehensive understanding 
of the protean rogue wave phenomenon is still far from 
being acheived 0, M . Indeed these waves not only appear 
in oceans but also in the atmosphere [9[, in optics [6[, in 
plasmas [lO|, in superfluids [11|, in Bose-Einstein con- 
densates [12] and capillary waves [i3|- Peculiar aspects 
and common features of the multifaceted manifestations 
of rogue waves in their different physical realms are a 
subject of intense scientific debate 0. New studies of 
rogue waves in any of these disciplines contribute to give 
a global view on a complex process that to a large extent 
remains unexplored [l^ . 

Nonlinear dynamics is one of the theoretical framework 
that has been successful in predicting the basic features of 
rogue waves [151, [i6| . A formal prototypical description of 
a single rogue wave is provided by the so-called Peregrine 
soliton, a solution of the focusing nonlinear Schrodinger 
equation (NLSE) [13, [HI which features a rational depen- 
dence on both space and time coordinates. Such solution 
describes the growing evolution of a small, localized per- 
turbation of a plane wave whose peak subsequently gets 
amplified by a maximal factor 3 over the background and 
eventually decays and vanishes. After decades of debate 
S 0] 7 the Peregrine soliton has been observed experimen- 
tally only very recently in fiber optics [l9| , in water-wave 
tanks [20|, and in plasmas [21|. Moreover the Peregrine 
soliton turns out to be just the first of an infinite hierar- 
chy of higher order rational solitons of the focusing NLSE 
with a progressively increasing peak amplitude. Again 
their amplitude over the background is expressed by the 
ratio of progressively higher degree polynomials and are 
therefore localized both in space and time [22|, l23| • These 
theoretical findings and experimental observations prove 



that the approach based on fundamental nonlinear mod- 
els, such as NLSE, may be fruitful and appropriate to 
rogue waves description. 

In a variety of physical contexts, several waves rather 
than a single one need to be considered, in order to ac- 
count for important resonant interaction processes. In 
these circumstances extreme waves should be described 
as solutions of coupled systems of equations rather than 
by the single-wave NLSE model. In this direction the 
investigation of solutions which are possible candidates 
as rogue waves has been recently extended to coupled 
NLSEs |23-[27|. This weak resonant interaction of two 
waves has been shown to cause wave behaviors which 
could not be detected by the single Peregrine soliton. 
However, in order to account for the appearance and dy- 
namics of extreme waves in strong resonant processes, the 
three-wave resonant interaction (TWRI) equations seem 
to be the fundamental and universal model. Indeed this 
system describes the propagation and mixing of waves 
in weakly nonlinear and dispersive media. Applications 
are found in fluid-dynamics (capillary-gravity waves, in- 
ternal gravity waves, surface and internal waves), in 
optics (parametric amplification, frequency conversion, 
stimulated Raman and Brillouin scattering), in plasmas 
(plasma instability, laser-plasma interactions, radio fre- 
quency heating), in acoustics and solid-state physics. 

In this Letter, we introduce, for the first time to 
our knowledge, a family of rational multi-component so- 
lutions of the TWRI equations which describe unique 
events, i.e. "amplitude peaks" which are distinctive of 
rogue waves as being much higher than the surrounding 
background and well isolated in both space and time. 
These solutions are expected to be crucial in forecast- 
ing and explaining extreme waves in a variety of multi- 
component resonant processes (f.i. oceanography, optics, 
plasma physics). 

TWRI equations and Rogue Waves. The system of 
three coupled partial differential equations we chose here 
to model the resonant interaction of three waves in 1 + 1 



dimensions reads as follows (in the notation of \2§\): 
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where each subscript variable stands for partial differ- 
entiation. En = En{z,t), n = 1,2,3, are complex am- 
plitudes, t is the evolution variable and z is a second 
independent variable. The coefficients Vn are the ve- 
locities of the three waves and we assume the ordering 
Vi > 1^2 > "^^ With no loss of generality, we assume 
Fs = by writing these equations ([!]) in the reference 
frame moving with the same velocity of Es. The signs of 
the coupling constants, with the minus sign only in the 
equation with the intermediate velocity V2, correspond to 
the so-called "soliton-exchange" case in the terminology 

of m. 

It should be pointed out that the meaning of the com- 
plex amplitude En, and of the coordinates t, z, depends 
on the particular applicative context (f.i., fluid-dynamics 
1301 . plasma physics [31], nonlinear optics [32], acoustics 

TWRI equations, like NLSE, possess rational solutions 
with the property of representing, in each of the three 
waves En^ "amplitude peaks" which are isolated in space 
and time. Similarly to the case of the NLSE, these so- 
lutions are local deformations of a non vanishing back- 
ground whose modulation instability is discussed in [3^ . 
Such solutions can be expressed as: 



El = 2q6i [1 ■ 
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i^2=2,5,[l + ^,^^^^|^J4^]e-(*+-^) (2b) 



E3 = 2iqS3[l 
where 
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(2c) 



= {-V3^i)/2,\0\ = l, 



Si = V{Vl-V2)/Vi,S2 = V(^l-l^2)/^2, 



S3 = {Vl - V2)/VV^2. 

^ = 71 + 726 + 73(^ - ^^*), 

^1 = 71 + 72(a + ^*) + 73(^ + 6>*a + iV3), 

^2 = 71 + 72(6 + 6>) + 73(7? + 6>^i), 

£,i=-2q{t^ipiz), r]= -^l-2iqp2Z, 

pi = e/Vi - r/F2, P2 = iM - 1/^2, 

zyi= 2/^2-1/^1, 

Z^2 = 1/^2 -2/^1. 



The above expressions depend on the velocities, Vi,V2, 
the real "frequency" parameter g, and the complex pa- 
rameters 71,72,73. Once the structural parameters (i.e. 
the characteristic velocities Vi and V2) are fixed, we are 
left with four independent parameters g, and 71,72,73. 





FIG. 1: Vector rogue waves envelope distributions |^i|, |^2| 
and l^sl of ([SI). Here, Vi = 1,^2 = 0.5,^ = 1,71 = 1,72 = 
1,73 = 0. 

However not all these parameters are essential, since 
some of them can be fixed without loosing generality by 
using appropriate symmetries of the TWRI equations ([T]) . 
The parameter q merely rescales the wave amplitudes, 
and the coordinates z and t. Thus one can set g = 1. 

Also the three remaining parameters 71 , 72 , 73 are not 
all essential as one (non vanishing) of them can be 
given the unit value. Moreover it can be shown that if 
72 = 73 = then the solution (|2]) represents plane wave 
backgrounds with no interest. Otherwise, if 73 = then 
the parameter 71 can be made to vanish by using trans- 
lation invar iance, while, by the same argument, one can 
set 72 = if 73 7^ 0, while 71 remains instead an essential 
parameter. Despite this simplification, we choose to keep 



all three parameters 71 , 72 , 73 and to play with them to 
better display a few aspects of the many properties of 
this family of solutions (J2]) . 
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FIG. 3: Vector rogue waves envelope distributions \Ei\, \E2\ 
and |^3| of 0. Here, Vi = 1,^2 = 0.5,^ = l,7i = 1,72 = 
0.1,73 = i. 



FIG. 2: Vector rogue waves envelope distributions \Ei\, \E2\ 
and l^sl of {|2)). Here, Vi = 1,^2 = 0.5, ^ = 1,71 = 2,72 = 

7,73 = 1.5 -\-i. 

In Fig. [U we first show the case 72 7^ 0, 73 = 0. The 
parameter 71 is so chosen as to put the peak at the origin 
of the (z, t) plane. As expected, the expression (|2]) de- 
scribes amplitude peaks which are localized in both z and 
t. Interestingly, each component \En\ looks like a rogue 
wave whose maximum height is twice the background 
intensity while its minimum is zero; its eye-shaped dis- 
tribution density shows one hump and two valleys. We 
also note that, as for the Peregrine soliton, the rational 
expression ([2]) is the ratio of two polynomials of second 
degree in the coordinates z^t. 

In the case 73 7^ 0, the expression (|2]) may describe 
amplitudes with multiple peaks localized in z and t. Fig- 
ure [2] shows two rogue waves with different structures in 
each one of the three components \En\. Figure [2] shows, in 
El component, a bright rogue wave with an eye-shaped 



distribution (a hump and two valley), together with a 
wave with a four-petaled distribution (two humps and 
two valleys around a center, and the center value is al- 
most equal to that of the background). The four-petaled 
wave in Ei component corresponds to eye-shaped rogue 
waves in E2 and Es components. By decreasing the value 
1 73/72 1, in each component these rogue waves separate. 
By increasing I73/72I, these rogue waves instead merge, 
giving birth to higher-amplitude vector rogue wave so- 
lutions (see Fig. [3j). The maximum value of the humps 
is more than 3 times the plane wave's background for 
some components, and the minimum value of the valleys 
is zero. 

Notice that effective energy exchanges take place be- 
tween waves ^1, E2 and Es in TWRI during their in- 
teraction. Figure [H reports a typical evolution of the 
effective energy /i, I2 and Is versus t. Effective energies 



.4^ 





FIG. 4: Effective energy evolution /i (blue line), I2 (green 
line) and 1 3 (red line) versus t, for the case reported in Fig. 
[I](left) and Fig. [3] (right). 



/i, I2 and /a are obtained according to the prescription: 



\E^ 



^nO| 



')dz, 



where E^q = lim^^oo^^n, ^ = 1^2, 3, is the plane wave 
background. The energy transfer between the waves can 
enhance the peak amplitude in some of the wave compo- 
nents. This wave behavior is completely different from 
what happens in coupled NLSEs, where energy exchanges 
are forbidden. 

Let us briefly discuss the experimental condition in 
nonlinear optics for the observation of TWRI rogue 
waves. In fact, nonlinear optics has been recently seen 
as a fertile, reproducible and safe ground to experimen- 
tally develop the knowledge of rogue waves [6|, 0, [ij, [l9| • 
One may consider a TWRI optical spatial non-collinear 
scheme with type II second-harmonic generation in a 3- 
cm long birefringent KTP crystal (f.i., see the experi- 
mental set-up of Ref. |32|]). Spatial diffraction- less 6-mm 
waist beams, mimicking quasi-plane waves, at 1064 nm 
(o-wave, and e-wave) and at 532 nm (e-wave) would lead 



to TWRI modulational instability evidence and rogue 
wave dynamics with peak field intensities of tens of 
MW/cm^. 

As a final remark, for the way of computing the expres- 
sion (j2j), we limit ourselves to notice that this is based 
on the Darboux technique applied to the Lax pair associ- 
ated to the TWRI equations. This method is well known 
and does not need to be detailed here to any extent. The 
relevant literature is rather vast and we refer to [35| for 
the formalism we have adopted and to [23, [36| for the 
basic arguments to follow for the construction of rational 
solutions. 

Conclusions. We have reported the explicit analytic 
expression of solutions of the equations describing the 
resonant interaction of three waves. These solutions 
have the important property of describing rogue wave 
events. Several articles have been recently devoted to 
rogue waves as rational solutions of multi-component sys - 
tems of coupled wave equations: VNLS equations [2J- 
27], Davey- Stewart son equation [STi] and coupled Hirota 
systems [38|. 

The present step in this direction, dealing with rogue 
wave solutions of the TWRI equations which represents 
a fundamental and universal model for the description of 
strong resonant interactions, seems to be a crucial stride 
to controlling and forecasting extreme- wave dynamics in 
multi-component wave systems, with a broad variety of 
applications. 
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